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The Einstein equations are completely integrated in the presence of two (incoming and outgoing) 
streams of null dust, under the assumptions of spherical symmetry and staticity. The solution is 
also written in double null and in radiation coordinates and it is reinterpreted as an anisotropic 
fluid. Interior matching with a static fluid and exterior matching with the Vaidya solution along 
null hypersurfaces is discussed. The connection with two-dimensional dilaton gravity is established. 
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I. INTRODUCTION 

Null dust represents the high frequency (geometrical optics) approximation to the unidirectional radial flow of 
unpolarizcd radiation. This is a reasonable approximation whenever the wavelength of the radiation is negligible 
compared to the curvature radius of the background. Various exact solutions of the Einstein field equations were 
found in the presence of pure null dust (for reviews see [jj], and more recently (|] and [[|). 
Oh' In some scenarios even gravitation behaves as null dust. Price has shown that a collapsing spheroid radiates 
away all of its initial characteristics excepting its mass, angular momentum and charge. (This result is known as 
the no hair theorem.) The escaping radiation then interacts with the curvature of the background being partially 
backscattered. Both the escaping and the backscattered radiation can be modeled by null dust |(| as the curvature 
radius of the background is larger than the wavelength of the radiation. 

Letelier has shown that the matter source composed of two null dust clouds can be interpreted as an anisotropic 
fluid 0], giving also the general solution for plane-symmetric anisotropic fluid with two null dust components. Later 
' Letelier and Wang [|| have discussed the collision of cylindrical null dust clouds. The collision of spherical null dust 
streams was discussed by Poisson and Israel Q. Their analysis yielded the phenomenon of mass inflation. 
0^ ' However, no exact solution in the presence of two colliding null dust streams with spherical symmetry was known 
until now. Recently Date tackled this problem under the assumption of staticity, integrating part of the Einstein 
equations. It is the purpose of the present paper to present the exact solution for the case of two colliding spherically 
symmetric null dust streams in equilibrium, in a completely integrated form. 

The plan of the paper is as follows. In Sec. II we derive the field equations and we integrate them. The emerging 
exact solution is written explicitly in suitable coordinates adapted to spherical symmetry and staticity. The metric 
Ch . in radiation coordinates and in double null coordinates is also given. We analyze the metric both analytically and by 
numeric plots. In Sec. Ill we present various possible interpretations of the solution, including the anisotropic fluid 
, picture, and dilatonic gravity. 

Finally Sec. IV contains the analysis of the interior and exterior matching conditions on junctions along timelike 
and null hypersurfaces, respectively. We employ the matching procedure of Barrabes and Israel po| . The interior 
junction fixes the parameters of the solution in terms of physical characteristics of the central star: the mass and 
energy density on the junction. The exterior matching with incoming and outgoing Vaidya solutions p]| leads to 
the conclusion that no distributional matter is present at the junction. This feature is in contrast with the exterior 
matching with the Schwarzschild solution H. 



II. SOLUTION OF THE EINSTEIN EQUATIONS 

The general form of a spherically symmetric, static metric in a spacetime with topology R 2 x S 2 is 

ds 2 = -h{r)dt 2 + /(r^dr 2 + r 2 dfl 2 . (2.1) 

Here t is time, r is the curvature coordinate (e.g. the radius of the sphere i=const with area Anr 2 ), dfl 2 = d9 2 + 
sin 2 (0)d(f) 2 is the square of the solid angle element. The functions /(r) and h(r) are positive valued. We introduce 
the local mass function m(r), related to the gravitational energy within the sphere of radius r ]l2| |: 

m = 1 - . (2.2) 

r 

The energy- momentum tensor in the region of the cross-flowing null dust is a superposition of the energy- momentum 
tensors of the incoming and outgoing components: 
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T ab = EfLf v a v b + u a u b\ 

8nr A 



(2.3) 



All energy conditions are satisfied for (3(r) > 0. The same linear mass density function (3 was chosen for both 
components as staticity requires no net flow in either of the null directions. The vector fields 



(2.4) 



are the tangents to the (future-oriented) outgoing and incoming null congruences, partially normalized such that 
v a u a — —1. Similarly as for the one-component null dust, here both null congruences are geodesies 0. 
After eliminating second time derivatives from the nontrivial Einstein equations we have the system: 
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(2.5) 



The solutions with /3=const all reduce to = 0, meaning vacuum. They are either the Schwarzschild or the flat 
solution, in accordance with the Birkhoff theorem. For (3 ^const one of the equations is a relation between the mass 
density 0(r) and the metric function h (r): 



h = 



A 

0' 



where A > is a constant. This relation can be deduced also from the energy-momentum conservation 
remaining equations do not contain the constant A: 
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Now we solve the system (2. 7), (2.. 



d\nr 

ncreas 
Following Date 



/3-/ + 1 
+ f-l 



(2.6) 
The 

(2.7) 
(2.8) 



Inserting (2.2) in (2.7) the m ass is fou nd to increase with the radius: dm/dr — (3/2 > . 

, we eliminate / from (2.7) by its expression taken from 



/ = 



l + (3 



1 - 



The resulting second order ordinary differential equation in 1//3 can be integrated, finding: 
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(2.9) 



(2.10) 



Here D is an integration constant. Inserting this expression in (|2.9| ), an algebraic relation between the metric compo- 
nents emerges: 



(1 + /3) 2 

V0 



(2.11) 



where C = exp((l + D)/2) > 0. 

Then we complete the integration of the system (2.7),(2.£). The key remark is that none of the equations in this 
system contains explicitly the independent variable lnr. Thus one can pass to the new independent variable 0, in 
terms of which an ordinary first order equation can be written: 



df = f(0 + f- 1) 
d(3 0(0 - f + 1) 



(2.12) 
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Introducing the new positive variables 

P= (2//?) 1/2 , 



L = 



1 + 



(VP) 1 ' 2 

the equation ( p. 12 ) takes the form of a first order linear (inhomogeneous) ordinary differential equation: 

dP 



dL 



2(1 - PL) 



The solution is found by integrating first the homogeneous equation, then varying the constant. It is 

r-L 



P = 2e- L2 <f> B (L) , $ B (L) = B + J e x2 dx>0, 



(2.13) 



(2.14) 



(2.15) 



where B is a third integration constant. The function &b(L) in ( 2.15 ) can be expressed either in terms of the error 
function or in terms of the Dawson function: 



&b(L)-B 



•erf(iL) = e L Dawson(L). 



(2.16) 



For properties of th e se tr anscen dent al functions see |l3[ ]. 

From Eqs. (2.11), ( 2T3| ) and (2.15) both the curvature coordinate r and the metric functions (5 and / are found as 
functions of the radial variable L: 



Cr(L) 



2L$ B (L) 



0(L) = -l + 2Le- L2 <Z> B {L) 
f(L) - »M 



\_ e V + 2L$ B (L)] • 



Then the mass function m = m (L) is obtained from (2.2): 



2Cm(L) = -e L2 + 2i$ B (L) - 2e" L2 $|(L). 



It is easy to check that both r and m are monotonously increasing functions of L: 

d(Cr) 

d(Cm) 



dL 



2<S> B (L) > 



dL 



<$> B {L) 



2L$ b (L) 



> 0. 



(2.17) 
(2.18) 

(2.19) 



(2.20) 

(2.21) 
(2.22) 



In the last relation the equality holds for r — 0. 

Now we have everything together to write the metric in terms of the new radial coordinate L: 



ds 2 



-Ae dt 2 



-e L2 +2L$ B {L) 



2e L 



2L$ B (£) 



dL 2 

-Q2 



2L4> B (L) 



2 dVL 2 



(2.23) 



There are three parameters in the solution, two of them restricted to be positive: A and C. Without loss of 
generality we can choose A — 1 by rescaling the time coordinate. The parameter C provides some distance scale. We 
comment on the third parameter in what follows. Both Cr > and the energy conditions [3 > imply 



B > X (L) , X (L) 



e 

~2L 



e x dx 



(2.24) 



valid for all admissible values of L. The equality B = x(Lq) holds for L = Lq corresponding to r = 0. As 
dx/dL — —e L /2L 2 < 0, the function x(L) is monotonously decreasing and the inequality ( 2.24 ) will be satisfied for 
any L > Lq. Thus B gives the lower boundary Lq for the range of the radial coordinate L. 

Next we plot numerically the functions /3(r), f(r) and m (r) for different values of the parameter B (Figs. 1-3). 
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FIG. 1. (a) The function (5 = /?(lnr) at parameter values B = C = 1. The metric is singular at r = and the space-time is 
not asymptotically flat, (b) Plot of f3 = (3(lnr, B) for C = 1 and < B < 7. 




FIG. 3. (a) The mass function m = m(lnr) at parameter values B = C = 1 takes negative values close to r = 0. (b) The 
mass function in the parameter region < B < 7. The value of r where the mass becomes negative depends on the parameter 
B. 



As we see from Fig. 3, the mass function vanishes at some radius f and takes negative values at r < f. We show 



here that such a positive f exists irrespective of the choice of the parameters B and C . By combining (2.17) and 
( p. 20 ) with the condition to = we find f = 2$^(L)/Cexp(L 2 ) > 0. The constant C has a simple scaling effect on 



the value of f. Figure 3(b) shows that in the domain (0, 7) the negative mass region can be extended by increasing 
B. We emphasize that the energy conditions are still satisfied in the negative mass regions. The interpretation of the 
negative mass function is not immediate and depends on the measuring procedure of the mass in this asymptotically 
non-flat space-time. 

At the end of this section we give the metric in double null and radiation coordinates. This requires an additional 
integration. We introduce a "tortoise coordinate" r* in the same manner it can be introduced in the Schwarzschild 
space-time: 



dr' 



[f(r')h(r')} 



1/2 



V2 
C 



2y<f> B (y) dy 



(2.25) 



The radial null geodesies are t = to + cr* , with c = 1 for incoming and c = — 1 for outgoing geodesies. 
Introducing the null coordinates x — t ± r* the metric takes the simple form: 

ds 2 = -h(x + ,x~)dx + dx~ + r 2 (x + , x~)dn 2 . 



(2.26) 



The radial null geodesies are now a: ± =c o nst. T he fun ctions h(x + , x~) — h(x + — x~) and r(x + , x~) — r{x + — x~) are 
contained in implicit form in Eqs. (2. 6), (2. 25), (2.17) and ( 2.18). 

Finally we cast our static, spherically symmetric solution of crossflowing null dust in either the incoming or outgoing 
radiation coordinates (v = t + cr* , r, 6, cf>). These coordinates are like the Eddington-Finkelstein coordinates for the 
Schwarzschild solution. If the solution was asymptotically flat, they would be the Bondi-Sachs coordinates. 



ds 1 



UP) 



-1/2 



dv 



1/2 



dv — 2cdr 



+ r 2 dfl 2 . 



(2.27) 



Here v = x for c = ±1. In these coordinates the radial null geodesies are given by one of the equations u=const and 



v = const + 2c 



Now it is evident that there is no apparent horizon: 



fir') 



1/2 



dr'. 



dr 1 // 



d{cv) 2 \(3 



1/2 



>o, 



(2.28) 



(2.29) 



thus no event horizon either. Thus the singularity in the origin r = is naked. This is very similar to the naked 
singularity of the negative mass Schwarzschild solution . 



III. INTERPRETATION OF THE SOLUTION 



A. 2D dilatonic model 

In this subsection we present a dilatonic model which in 4D has the interpretation of a spherically symmetric 
gravitational field in the presence of two crossflowing null dust streams. This dilatonic model emerges from the action 



5' 



d 2 Xy 



-99 



1 



(3.1) 



The first term i s th e Einstein-Hilbert action reduced by spherical symmetry (a surface term was dropped). The 
second term in (3.1) represents a 2D massless scalar field in minimal coupling. The minus sign assures that all 4D 



energy conditions are satisfied, p = r 2 is the dilaton and ip is the sc alar field. The conformal flatness of the 2x2 



metric g a /3 — hn a p is manifest in the corresponding 4D line element ( 2.26| ). V and R [g] are the covariant derivative 
and curvature scalar, respectively associated with the metric g. Although the equations emerging from this model 
are quite similar to the corresponding equations in the Callan-Giddings-Harvey-Strominger (CGHS) model E3], they 
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could not be exactly integrated in general. For our purposes we take only the equation emerging from variation of 
the scalar field, in double null coordinates: 



*»,+-= 

Here commas denote derivatives. The equation has the D'Alembert solution 

ip = <p+ (x+) + <p~ (x~) , (3.2) 

showing that the scalar field behaves like our matter source of crossflowing null dust. The 4D interpretation of the 
particular solution with either leftmoving or rightmoving matter is the Vaidya solution JTlJ] . Recently Mikovic has 
given a solution |l5| for the case where both components are present, in the form of a perturbative series in powers 
of the outgoing energy- mom entum component. 



Our static solution ( 2.26 ) represents the first explicit exact solution for this model, when none of the null dust 



components are neglected. 



B. Anisotropic fluid 



Following Letelier J7[ we reinterpret the energy- momentum tensor (2^3) as describing an anisotropic fluid with a 
pressure component equaling its energy density: 

Tab = p(U a U b + XaXb) (3.3) 
~U a U a = X aX a = 1 , U aX a = . (3.4) 



A straightforward comparison with (J2 .3|) gives 



" = ^7" (3 - 5) 
i . . i fsv 



-*>-^(£)'- (3 - 7) 

Thus the solution represents an anisotropic fluid at rest with the energy density p and radial pressure p = p. No 
tangential pressures in the spheres r=const are present. The fluid is isotropic only about a single point, the origin. 

C. Radiation atmosphere 

The solution can be interpreted as the outer region of a radiating star, receiving radiation from the s urrounding 



region either. If equilibrium is achieved between the two components, we have the static solution (2.23). This was 
the initial interpretation proposed by Date [S . 

We write the energy-momentum tensor of the static solution in the double null coordinate system (x + , x~, 6, <p). 
Inserting the null covectors 

/ fu \ I Fu \ 

(3.8) 





in the covariant form of (|2.3|) and taking account of (2.6) the energy- momentum tensor becomes 



T ab dx a dx b = — i-^r (dx+dx + + dx~dx~) , (3.9) 

a superposition of two cross-flowing null dust streams with equal and constant mass density functions. However one 
can freely rescale the null vectors v a and u a to have arbitrary mass density functions either. After such a rescaling 
the mass density function of the incoming null dust depends only on the outgoing coordinate x + and viceversa, a 
property pertinent to the mass functions of the Vaidya solution , characterized by the energy- momentum tensor 



6 



T ab dx a dx b = ° dM Jp dVdV. (3.10) 
4irr z dv 

M(V) is the mass function and the coordinate V is outgoing for incoming radiation and incoming for outgoing 
radiation. 

In the next section we will study the interior junction with a static star, and the exterior junction with incoming 
and outgoing Vaidya solutions. 

IV. JUNCTION CONDITIONS 

A. Matching with interior spherically symmetric static solutions 

We discuss the junction with an interior solution with accent on the anisotropic fluid interpretation given previously. 
A similar treatment was given in Both our analysis and the one in || reveals that the junction with a static 
interior matter can be done without a regularizing thin shell. Our treatment is more general, however. We formulate 
the junction conditions for two generic static spherically symmetric space-times, following the standard Darmois-Israel 
junction procedure [[l6| , fL7| and we establish a constraint on the matter pressures implied by the matching conditions. 
An other improvement over II is due to the fact that we dispose of the exact solution ( 2.23Q , thus the explicit 



computation of the the matching conditions with an arbitrary particular interior becomes possible. 



An orthonormal basis is given by the vectors U and \ defined by Eqs. ( [3.6] ) and (3.7) together with the spacelike 
vectors 



ld_ _ 1 3 

r 39 ' 4 r sin 9 36 



(4.1) 



Any spherically symmetric static energy-momentum tensor has the form 

T ab = P U a U b + PlX a X b + P2 {E a 3 E b 3 + E%E b A ) . (4.2) 



By inserting pi = cul^ixr and p in the form (3.5) in the Einstein equations for the metric (2.1) we find 



dr 

rf^ = h( ai -f + l) (4.3) 

2rf ^dV = ~ a ' + ai{ ~ fj + / - 1) + w - !) + 4 "2/- 

The induced metric of the surface r=const is 

ds| = -hdt 2 +r 2 dn 2 . (4.4) 

Without loss of generality we can choose the time coordinates in both static space-times such that they are continuous 
on the junction. Then the continuity of the first fundamental form requires the metric function h to be continuous. 
The extrinsic curvature of the junction surface r = const is defined as 

K ab = (S c a - X c Xa) {Si - X d Xb) V cX d . (4.5) 

The nonvanishing components are: 

Jf tt = * (- ai + /-l), Kl=Kl = ^L. (4.6) 
2rv/ r 



In the above expressions the derivatives were eliminated by use of Eq. (4.3). 

Continuity of the extrinsic curvature across the junction hypersurface r = r% is achieved, provided that the metric 
and the function ct\ (thus also p±) are continuous. We have proved the following result: 

Any two spherically symmetric static solutions can be matched along hypersurfaces r = const provided the radial 
pressures are continuous. This is similar to the theorem given by Fayos, Jaen, Llanta and Senovilla [|18[ for the 
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matching of the Vaidya solution with a generic spherically symmetric solution along timelike hypersurfaces. A generic 
discussion on matching spherically symmetric space-times along thin spherical timelike shells can be found in fl9| . 

For the double nu ll d ust solution a-y — (3. In consequence the interior fluid should have a radial pressure equal to 
the energy density (3.5) of the double null dust solution on the junction. However, no conditions on the pressures 
tangent to the spher es em er ge. 

We see from Eqs. ( |2.18 ), ( 2.20 ) and ( j3~5| ) that the integration constants B and C appear in the radial pressure and 
mass function of the static double null dust solution. Continuity of these functions on the junction fixes the value of 
the constants, once the interior solution is chosen. For a realistic star, the mass Mi should be positive. This implies 
a lower boundary for the possible values of n, as follows from Fig. 3(a). 



Let us illustrate the junction with the interior Schwarzschild solution [EOl ds 2 



(a - by/1 -r 2 /i? 2 ) dt 2 + 



'l-r 2 /R 2 



dr 2 



2 df2 2 , with the energy density p$ = 3/i? 2 =const and pressure pg given by 



3b Jl 



8nps = 



R 2 la-bJl 



(4.7) 



Several relations among the Schwarzschild parameters a, b, R, 
emerge from the junction conditions: 



radius of the star r± and the parameters B and C 



2mi 



R 2 ' 



2a 



R/_ 

r 2 ' 



'Hi , 



2b\ 1 



' l 
R 2 



R/_ 

r? ' 



(4.8) 



We have denoted by fix and mi the values of the functions (3 and m at the junction r = r±. The first two relations ( |4.8j ) 
determine the constant B and the value of the radial coordinate L\ at the junction in terms of r\jR and a, when Eqs. 
( |2.17 ), (2.18) and (2.20) are inserted. Eliminating Pi from the last two r elati ons of d4.8| ), a constraint b = b(a,ri/R) 
on the possible values of characteristics of the interior emerges. Finally ( [ml ) implies C = C(a, n/R, r x ). 

In the light of the above relations we see that after choosing some value for C (a scale), the constant B is determined 
exclusively by the radius and density (or mass) of the star. 



B. Matching with Vaidya solutions 

In this subsection we study the junctions with the incoming and outgoing Vaidya solutions, which are at the 
exterior of the static double null dust solution. There are only three points (in fact spheres) in common with exterior 
Schwarzschild regions (Fig. 4), and the matching can be done without introducing regularizing thin shells. 

The high-frequency approximation to the uni direct ional radial flow of spherically symmetric unpolarized radiation, 
characterized by the energy- momentum tensor ( [3.10 ) is represented by the Vaidya solution JO]: 



ds 2 =-dV 



The radial null geodesies are the lines V=const and the curves given by 



r 2 dfl z . (4.9) 



* = c - ( i - ™m) . (4.10) 

dV 2 V r J y ' 

We would like to match our static solution with incoming and outgoing Vaidya solutions along the outgoing, 
respectively incoming radial null geod esies, e.g. along lines of constant in comi ng, respectively outgoing coordinate 



(Fig 4.). These are given by Eq. ( 4.10 ) in the Vaidya solution and by Eq. ( p.28| ) in the static solution. 

A convenient formalism for matching solutions along null surfaces, which does not require coordinates that match 
continuously on the shell, was developed by Barrabes and Israel pToj ] . Their discussion on the particular case of 
spherical symmetry requires the metric in both space-times written in the form: 

ds 2 = -e^dz [fe^dz - 2dr] + r 2 dtt 2 . (4.11) 

Here / and ip depend on both r and z. The null coordinate z is always outgoing but it can be either increasing or 
decreasing with time. Then the junction is done along the null hypersurfaces z=const. 
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To apply this formalism it wou ld be desirable to express the Vaidya solution in the other set of radiation coordinates 
[U(V,r),r,8,cj)} in which Eq. flyg ) takes the form U =const. However, this is possible only when double null 
coordinates are found. It was demonstrated in plj that double null coordinates exist when the mass is a linear or 
exponential function of the advanced or retarded time V. We argue in what follows that the mass function of that 
particular Vaidya solution, which is matching continuously to the static solution has a more complicated dependence. 
A further inconvenience is that there is no obvious choice for the intrinsic coordinates in terms of the space-time 
coordinates in the null junction hypersurface. For these reasons we proceed as follows. First we find coordinates 
that match continuously on the junction and in terms of which the metric is continuous. Then we employ the 
Barrabes-Israel junction formalism to find the distributional stress-energy tensor on the junction. 




FIG. 4. The intersection of two crossflowing null dust streams with mass functions M(V) is the static solution characterized 
by the mass function m(r). Here V = X stands for c = ±1 and X are the advanced (retarded) time in the Vaidya solution. 
The 2-dust solution touches in three points exterior Schwarzschild regions with masses Mi and M2. At the interior junction 
there is a static fluid representing a star with mass Mi . The mass functions M(V) of the Vaidya regions change monotonously 
from Mi to M 2 . 



As a radial variable of both metrics we choose L by extending the expression (2.17) of r = r(L) to the Vaidya 
regions too. An appropriate null coordinate z is defined by the values of L = Lj on the junction. In these coordinates 
the junction hypersurfaces are characterized by the null geodesic equations z = L in both space-times. We proceed 
in deriving the expressions of the coordinates v and V and of the mass function M in terms of the coordinate z. 

Identifying the corresponding part of Eq. ( |3.9[ ) with Eq. (3.10 ), we find the relation between the null coordinates 
V and v, valid on the junction : 



Then we extend this relation over b oth the static and the Vaidya regions. 

The null geodesic equation ( 2.28 ) of the static solution, evaluated on the junction (where z = L) gives 

dz _ c dM 
dV ~ 2r 2 dz 



(4.12) 



(4.13) 



Inserting Eq. ( 4.13 ) in the null geodesic equation ( 4.10 ) of the Vaidya solution we find its mass function in terms of 
the new null coordinate z: 



2CM(z) 



2z^ B {z)~2e~ z $|(z). 



(4.14) 



Thus the mass function is continuous at the junction. The geodesic equation (4.13) together with (4.14) gives the 
relation between the null coordinates V and z: 
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2c 

V (z) — const + — 



2y<$> B {y) 



(4.15) 



The relations ( 4.14 ) and ( 4.15 ) contain in implicit form the dependence of the mass function M of the Vaidya solution 
on the radiation coordinate V. Despite the complicated dependence M (V) it is straightforward to check that M 
satisfies the required monotonicity condition: 



dM 
l d(cV) 



„2z 2 



> 



Finally Eqs. ( J4.12|) and (|4.15|) give v = v(z) 



v (z) — const + 



2V2c 
C 



2y^ B (y) dy 



Now we express both metrics in terms of the coordinates (z,r(L),6,<p). They take the form (4.11) with 



/=1- 
fv = l~ 



2m (L) 
2M (z) 

T(ZT 



,1> - 



,ll>V — 



C<P B {L) 

2 e z2 



-e z + 2z$ B {z) 
-e z2 +2z<f> B (z) 



(4.16) 

(4.17) 

(4.18) 
(4.19) 



Here the in dex V refers to the Vaidya solution and the expressions r (L) , to (L) and M (z) are given by Eqs. ( 2.17 ), 
(2.20) and (4.14), respectively. We have completed the task of writing both metrics in coordinates which are continuous 
on the junction and in terms of which the metric is continuous. 

It is immediate to check the continuity of the induced metric given by ds^ — r 2 dfl 2 . The other junction condition 
is a somewhat subtle issue as the conventional extrinsic curvature tensor for null hypcrsurfaces carries no transversal 
information. 



We define a pseudo-orthonormal basis |22| (n, N, £3, £4), where £3 and £4 are given in Eq. (4.1) and 



d f d 



dz 



2 dr 



N 



d_ 

dr 



(4.20) 



The vector n a is orthogonal (and also tangent) to the hypersurfaces fe^dz — 2dr = 0, along which the two space- 
times are glued together. The vector N a is the other radial null vector, transversal to these surfaces. They are 
related to the previously introduced vectors v a and 11° as follows. For the junction with the incoming Vaidya region 
v a = (2/ Z) 1 ^ 2 n a and u a = (f/2) 1 ^ 2 N a , while for the junction with the outgoing Vaidya region u a = — (2//) 
and v a = - {f/2) 1/2 N a . 

The projector to these null hypersurfaces with tangent space spanned by £3, £4 and n is 



1/2, 



KEt 



n a N b = 51 + N a n b . (4.21) 
Following we define the transverse or oblique extrinsic curvature tensor with the aid of the transverse vector N: 

K ab = P^V c N d . (4.22) 
The only nonvanishing components of the transverse extrinsic curvature tensor^] are 



C 



Kzz = ~^d r (/e*) = - 



, dL 
d(Cr) 



(4.23) 



Thus when matching any two metrics of the form (4.11) on the null hypersurfaces fe^dz — 2dr = 0, the jump in the 
extrinsic curvature is given by the jump in d r (/e 1 ^), provided the metric is continuous across the junction. 



x The components of the extrinsic curvature tensor defined in Jl^] are found from K a t by contracting with the three basis 
vectors E3, E4 and n tangent to the hypersurface. 
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A straightforward computation employing (4.18) and (4.19) gives on the junction hypersurfaces of the static double 
null dust region with the Vaidya regions 



Ol (/e ) = d L (f v e^) 



T 2 


T 2 


+ 2L^ B {L 3 )-e- L 




C 


'-e L *+2L ] <S> B (L j ) 





(4.24) 



In conclusion the extrinsic curvature is also continuous on the junction. There is no need for a thin regularizing 
shell separating the two domains of the space-time, in contrast to the exterior junction proposed in 0. 



V. CONCLUDING REMARKS 



We have integrated the Einstein equations in the presence of crossflowing null dust under the assumptions of 
spherical symmetry and staticity and analyzed various aspects related to the properties of the emerging exact solution 
(2.23). The solution has dilatonic gravity connections and it can be reinterpreted as an anisotropic fluid with radial 
pressure equal to its energy density and no pressures along the spheres r=const. This can be a radiation atmosphere 
for a star with its radial pressure equal to the energy density of the athmosphere on its surface. No constraint on the 
pressures along the spherical junction surface was found. On the exterior, the study of the matching conditions with 
the Vaidya solution revealed no thin shells on the junction. 

As a byproduct, we have der ived general conditions for the junction of two spherically symmetric solutions. Matching 
of two static space-times (2.1) along r=const hypersurfaces is assured by the continuity of the metric functions and 
of the radial pressure. Matching of generic spherically symmetric space-times ( [4.11 ) along the null hypersurfaces 
fe^dz — 2dr = is possible whenever the metric and d r (/e^) are continuous. 

The negativity of the mass function in some neighbourhood of the r — singularity raises the possibility of matching 
this solution to a negative mass core. This may be difficult due to the nontrivial topology of the known negative mass 
black hole solutions p3[ |. Despite the lack of experimental evidence for negative mass objects, presumably of quantum 
origin |]] , their microlensing effect |24|] on radiation from Active Galactic Nuclei was shown to produce features similar 
to some observed Gamma Ray Bursts p5[ . 

Equally interesting would be to introduce dynamics in the picture by an interior matching with a collapsing star. 
We defer this topic to a forthcoming study. 

An intriguing open question remains whether exact solutions describing the collision of spherically symmetric null 
dust streams, which have not reached equilibrium, can be found. 

Note added. After the submission of this paper a relevan t wor k in the subject was published by Kramer 
solution presented there is the particular case of the metric ( 2.23| ) with the parameter values B = and C = 
metric functions h and /, radial variable L and transcendental function — e L ~ +2L$ B correspond to e 2v ', e~ 2A , 
and y/ e(l + 2£) J, respectively of this paper, when the parameters valu es B = and C = ^fe are chosen. 




Keeping 

the parameters arbitrary enabled us in Sec. IV A to match the solution ( 2.23] ) with an interior Schwarzschild solution 
with arbitrary mass and radius. In contrast with our analysis relying on the Barrabes-Israel matching procedure, 
in p6| the junction with the Vaidya space-time was discussed by imposing the continuity of the four-metric on the 
junction. 
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